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13.  ABSTRACT 


The  research  program  in  Historical  Climatology  has  the  general  pur¬ 
pose  to  show  feasibility  of  determining  global  climate  in  the  past,  by 
determining  past  trends  in  surface  temperature  of  the  earth.  The  method 
being  studied  involves  measuring  past  temperatures  using  isotope  "thermom¬ 
eters"  in  bio-organic  data  bases.  For  example,  isotope  ratios  in  fee 
rings  belof  •foriiTO  are  expected  to  change  if  the  temperature  .changes ; 
thus  tree  rings  formed  long  ago  may  Cj&ibit  different  ratios  than  those 
being  fernod  today.  Lower  limits  on  the  temperature  coefficients  have 
been  computed,  and  compared  favorably  with  experimental  coefficients. 

New  measurements  of  these  coefficients  have  been  undertaken  in  water  plants. 
Feasibility  of  measurement  of  the  isotope  ratio  by  neutron  activation  for 
oxygen  has  been  demonstrated  and  feasibility  of  measurement  of  D/H,  C  /Zr 
and  0l8/0l6  has  been  demonstrated  by  elastic  backscatter  of  low  energy 
protons.  The  variations  of  isotope  reties  in  two  trees  have  been  measured 
and  the  correlation  with  historic  weather  records  is  being  studied*  The 
principle  and  application  of  multiple  thermometry  is  being  patented. 
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ANNUAL  TECHN7.CAL  REPORT 
PERIOD  ENDING  SEPTEMBER  30,  1972 
CONTRACT  NO.  F44620-72-C-0025 


The  research  program  in  Historical  Climatology  has  the  general  pur¬ 
pose  to  show  feasibility  of  determining  the  global  climate  in  the  past, 
determining  past  trends  in  surface  temperature  of  the  earth.  The  method 
being  studied  involves  measuring  past  temperatures  using  "isotope  ther¬ 
mometers,"  which  may  be  understood  in  the  following  way.  When  bio-organic 
material  is  formed,  such  as  cellulose,  in  trees,  the  isotope  ratios  of 
elements  in  the  bio-organic  material  may  differ  from  those  in  the  inor¬ 
ganic  material  from  which  the  living  organism  was  formed  and,  if  so,  the 
difference  depends  on  the  local  temperature.  For  example,  the  stable 
isotope  ratios  of  C^^/C^^  and  D/1I  in  cellulose  may  differ  as  a 

function  of  the  radius  of  the  ring  or  age  of  the  tree,  indicating  that 
from  time  to  time  in  the  life  of  the  tree,  the  climate  varied. 

As  a  first  step  in  initiating  this  research,  computation  was  made 
of  the  magnitude  of  effect  to  be  expected;  that  is,  the  temperature 
coefficients  of  isotope  ratios  of  oxygen,  hydrogen,  and  carbon  were  com¬ 
puted  in  cellulose  assuming  that  it  is  manufactured  at  equilibrium  with 
water,  carbon  dioxide,  and  oxygen  in  the  surroundings.  If,  instead, 
manufacture  does  not  occur  in  an  equilibrium  process,  then  the  actual 
temperature  coefficients  should  be  larger  than  those  computed  for  equi¬ 
librium.  Thus,  from  this  calculation,  we  should  obtain  lower  limits  of 
the  effects  to  be  expected.  The  report  of  these  calculations,  entitled 
"Multiple  Thermometry  in  Paleo  Climate  and  Historic  Climrte,"  by  Leona 
Libby,  January  1972,  has  been  published  in  the  Journal  of  Geophysical 
Research  77,  4310-17,  (1972). 
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The  temperature  coefficient  for  C  /C  in  cellulose  was  computed 
to  be  0.36  parts  per  thousand  per  degree  C  for  equilibrium  formation.  It 
is  in  agreement  with  an  experimental  value  measured  in  the  laboratory  for 
water  plants  by  Degans,  et  al.,  "Metabolic  Fractionation  of  Carbon  Isotopes 
in  Marine  Plankton,"  Deep  Sea  Research,  15,  1-9  (1968). 

The  temperature  coefficient  for  D/H  was  computed  to  be  2.4  parts  per 
thousand  per  degree  C.  Shortly  after  completion  of  the  calculation,  data 
from  which  an  experimental  value  for  D/H  in  peat  could  be  estimated  was 
published  by  W.  E.  Schiegl,  "Deuterium  Content  of  Peat  as  a  Paleo-Climate 
Recorder,"  Science  17_5 ,  512-513  (1972),  amounting  to  3.0  parts  per  thousand 
per  degree  C,  compatible  with  the  calculated  value. 
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With  the  help  of  Louis  Pandolfi,  a  graduate  student  in  the  Chemistry 
Department  of  UCLA,  an  experiment  to  determine  the  temperature  coefficients 
in  water  plants  has  been  undertaken.  In  an  outdoor  pond  of  800  liters 
volume,  water  plants  were  grown  and  the  nev;  growth  was  harvested  at  various 
constant  temperatures.  Dried  of  external  moisture,  the  plant  material 
was  freeze-dried  in  a  vacuum  system,  and  analyzed  for  isotope  ratios  of 
carbon,  hydrogen  and  oxygen  by  mass  spectrographic  methods.  The  measure¬ 
ments  are  being  analyzed  and  a  report  will  be  prepared  for  publication 
when  the  analysis  is  complete,  expected  to  occur  in  the  second  year  of 
this  research. 

We  have  analyzed  the  tree  rings  of  an  American  Bristlecone  Pine  and 
of  a  European  Oak  in  a  similar  way  for  stable  isotope  ratios,  and  we  are 
presently  undertaking  correlation  of  the  observed  variations  with  historic 
weather  records  at  stations  near  where  these  trees  grew.  A  report  on  this 
work  will  be  prepared  for  publication  when  it  is  complete. 

A  second  part  of  this  research  interests  itself  in  the  possibility 
of  determining  isotope  ratios  by  neutron  activation,  as  an  alternative 
to  mass  spectroscopy.  In  the  latter  method,  before  the  isotope  ratio  can 
be  measured,  many  chemical  reactions  are  involved  in  any  one  of  which 
isotope  fractionation  might  occur.  Thus,  using  an  alternate  method  for 
determination  of  the  ratio  would  be  important  as  a  check  on  the  method 
of  mass  spectroscope. 

To  this  end,  determination  of  the  ratio  C)18/ol&  to  the  required 
accuracy  has  been  demonstrated  to  be  feasible  by  measuring  the  activities: 

018 (n , y)018 (200KeVy) 

016(n,p)Nl6(6.1MeVy) 

induced  by  neutrons,  the  first  by  slow  neutrons  and  the  second  by  neutrons 
of  more  than  10  MeV.  Activation  should  be  made  in  a  reactor  with  a  suf¬ 
ficiently  high  flux  of  fast  neutrons,  such  as  the  Gulf  General  Atomic 
TRIGA  reactor.  This  method  is  non-destructive  of  the  sample.  Likewise, 
feasibility  of  measurement  of  the  stable  isotope  ratios  of  carbon,  hydrogen, 
and  oxygen  to  sufficient  accuracy  by  elastic  back-scattering  of  low  energy 
protons  has  been  demonstrated.  If  this  latter  method  is  to  be  non-destruc¬ 
tive,  then  it  requires  30  keV  protons  from  a  tandem  Van  der  Graaf  so  that 
reasonably  thick  targets  can  be  used. 
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The  ratio  0  /0  has  been  measured  to  6  parts  in  10,000  by  this 

method  in  about  40  minutes.  It  is  expected  that  soon  the  method  will  be 
reduced  to  a  routine  measurement  in  shorter  times  and  with,  perhaps,  some¬ 
what  better  accuracy.  A  report  is  being  prepared  for  publication. 
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A  patent  is  being  applied  for  concerning  the  principle  of  multiple 
thermometry  of  stable  isotopes  as  a  research  tool  to  study  past  climates, 
and  to  identify  and  correlate  sedimentary  layers  on  a  world  wide  basis. 
More  particularly,  we  intend  to  patent  the  use  of  this  principle  for  oil 
prospecting  and  prospecting  of  uranium  and  radium,  and  for  mapping  defor¬ 
mations  of  sedimentary  layers  caused  by  ice  loading  and  mountain  building. 
The  patent  search  has  been  completed  and  the  patent  instrument  is  being 
prepared. 
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Appendix  I : 

Multiple  Thermometry  in  Paleo  Climate  and  Historic  Climate,  by 
L.  M.  Libby,  J.  Geo.  Phy.  Res.  77,  No.  23,  4310-4317  (1972). 


Appendix  II: 
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Feasibility  of  Measurement  of  0  /0  Ratio  by  Activation  Using 
a  Fast  Neutron  Reactor,  by  S.  F.  Banker t,  S.  D.  Bloom,  L.  M.  Libby 
and  H.  R.  Lukens. 


Appendix  III: 

Determination  of  Isotopic  Abundances  of  C,  N,  and  O  in  Organic 
Samples  Using  Nuclear  Techniques  by  K.  W.  Jones  and  H.  W.  Kraner. 
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ABSTRACT 


The  principle  of  isotope  thermometry  is  that  isotope  ratios  are 
temperature  dependent.  It  is  proposed  that  isotope  ratios  in  bio- 
organic  material  may  have  independent  temperature  coefficients,  so 
that  by  measuring  ratios  for  several  elements,  e.g.  hydrogen,  carbon 
and  oxygen,  it  may  be  possible  to  show  that  a  temperature  change 
occurred.  Assuming  that  cellulose  is  formed  almost  at  equilibrium, 
temperature  coefficients  for  C,  H  and  0  are  computed  and  found  to  be 
of  measurable  size.  The  coefficients  computed  for  C  /C  and  D/H 
agree  with  measured  values.  Thus  indications  are  that  temperature 
changes  in  past  climates  may  be  measured  in  old  tree  rings,  using 
this  principle  of  multiple  thermometers. 


I.  INTRODUCTION 


The  principle  of  isotope  thermometry  is  that  the  isotopic  ratios 
are  temperature  dependent.  There  are  two  kinds  of  temperature  dependence, 
kinetic  and  equilibrium.  The  theory  of  kinetic  isotope  separation  is 
well  established.  For  example,  in  the  gaseous  diffusion  process  for 
separation  of  the  uranium  isotopes,  separation  is  caused  by  a  difference 
in  the  rates  with  which  the  two  species  pass  through  holes  in  a  barrier. 

In  case  of  equilibrium,  separation  is  caused  by  a  difference  of 

erences,  namely  the  difference  between  separation  by  the  forward 
rate  and  separation  by  the  reverse  rate,  the  two  rates  being  equal. 
Distillation  of  liquid  in  a  closed  container  is  an  example  of  separation 
at  equilibrium,  while  vaporization  into  a  vacuum  is  kinetic. 

In  distillation  of  a  liquid  containing  two  or  more  elements  in 
the  molecule,  the  corresponding  isotope  thermometers  are  functionally 
related  because  the  corresponding  vapor  pressure  ratios,  which  determine 
the  separations  on  evaporation  and  condensation,  depend  on  the  same 
force  constants  tying  them  together. 

In  bio-organic  material  however,  there  may  be  at  least  three 
useful  isotope  thermometers,  namely  D/H,  C13/C12,  and  018/016.  The 
basic  reaction  is  photosynthesis  with  production  of.  cellulose,  here 
written  schematically  as  a  basic  module,  H  -  C  -  OH,  or  CH20,  according 

^2  ®2*  reaction  all  three  elements  may 

yield  isotope  thermometers,  and  they  may  be  independent  of  each  other. 
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That  is,  the  ratio  0  / 0  is  no  longer  functionally  related  to  D/H 

because  diatomic  oxygen  is  being  evolved  and  so  rotational  and  vibra¬ 
tional  energies  of  the  0-0  bond  are  involved  and  because  in  cellulose 


-2- 


oxygen  is  bound  to  carbon  as  well  as  hydrogen.  For  the  sane  reason, 

the  temperature  coefficient  of  018/016  and  C13/C12  may  be  not  functionally 
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related.  And  finally,  the  ratio  C  /C  should  be  independent 

of  D/H  because  carbon  and  hydrogen  are  initially  in  different  molecules, 

with  differing  vibrational  frequencies. 


In  case  variations  are  measured  for  only  one  isotope  ratio  the 
relation  to' paleo  temperature  is  not  proven..  But  if  ratios  are 
measured  for  two,  three,  or  more  independent  isotopes  in  the  same  organic 
material,  it  is  possible  that  the  relation  to  temperature  can  be 
firmly  established  and  quantified,  and  this  is  the  principle  point 
of  the  present  discussion. 

Over  the  small  temperature  interval  (  0°C  <  t  <  +30°C)  within 
which  climate  could  have  varied  without  killing  terrestrial  life  as 
we  know  it,  the  dependences  of  isotope  ration  on  temperature  may  be 
assumed  to  be  linear  for  bio-organic  material.  That  is,  for  isotope 
ratios  fiJ,  where  j  =  2,  13,  15,  18,  34,  41,  .  .  .referring  to  D,  C13, 

N*3,  0*8,  S3\  K  .  .the  temperature  dependence  may  be  written  as 


^  s  t  4*  b^  , 


(a*5  and  bJ  constants) 


where  is  the  temperature  coefficient  and  t  is  temperature.  Then,  a 
change  At  is  related  to  a  change  AS3  by,  (for  a  given  compound), 


lt2  13  15  18 

At  =  =  AS.  =  ALL.  =  M_  = 


a. 


a13  a15  a18 
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The  multiple  over-determination  of  temperature  by  using  a  set  of  many 
thermometers  may  make  it  possible  to  show  with  some  level  of 
confidence  that  a  temperature  change  did  occur.  It  is  also  possible 
that  new  affects,  not  caused  by  temperature,  will  be  demonstrated  by 
the  method  of  multiple  thermometry,  if  some  of  the  above  ratios,  but 
not  all  of  them,  show  equality.  Probably  the  specific  relations  to 
temperature  will  have  to  be  calibrated  at  each  geographical  position 
to  take  into  account  local  variations  in  water  and  C02. 

Finally,  bio-organic  data  bases  can  be  radiocarbon  dated,  so  that 
if  evidence  for  temperature  changes  is  stored  in  them,  the  dates  when 
the  changes  occurred  can  be  evaluated. 


II. 


EXPERIMENTAL  EVIDENCE  FOR  TEMPERATURE  DEPENDENCE  OF 


ISOTOPE  RATIOS  IN  BIO- ORGANIC  MATERIAL 


Although  animal  bio-organic  material  is  not  directly  derived  from 
photo-synthesis,  nevertheless  animals  are  largely  what  they  eat,  so 


that  bio-organic  material  from  both  plants  and  animals  may  be  useful 


in  isotope  thermometry.  Some  of  the  evidence^^  on  isotope  thermometers 


in  bio-organic  material  is  as  follows. 
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The  C  /C  ratio  in  combustible  organic  matter  in  cores  from 


sea  bottoms  shows  variations  of  as  much  as  6  parts  per  thousand  with 


depth.  These  have  been  measured  and  attributed  to  temperature  variation 


in  the  sea  surface  with  time,  by  Rogers  and  Koons.^  Also,  in  a  study  of 


marine  plankton,  by  Sackett  et  al,^^  the  concentration  was  found 


depleted  by  6  per  mil  where  surface  waters  are  near  0°C  relative  to 
samples  collected  where  surface  temperatures  are  about  25°C.  These 


-TlVv 


temperature  coefficient  of  C  /C  in  bio-organT 


parts  per  thousand  per  degree  C  in  the  temperature  range  0  -  25°C. 

(A) 


In  further  work  by  Degens  et  al,  the  temperature  coefficient 
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of  C  /C  in  marine  phytoplankton  grown  in  laboratory  cultures  was 


measured  as  —0.35  parts  per  thousand  per  degree  C  between  10  -  30°C. 


While  shell  carbonate  is  not,  strictly  speaking,  bio-organic 

material,  nevertheless,  for  completeness  we  note  here  that  variations 
18  16 

in  the  0  /0  ratio  in  foraminifera  shells  from  Caribbean  Sea  cores 

have  been  measured  by  several  workers  and  interpreted  in  terms  of  paleo- 

temperature  variations. ^  This  in  fact  was  the  first  isotope 

thermometer,  as  developed  by  H.  C.  Urey.^^ 

It  may  be  that  other  isotope  ratios  in  bio-organic  material  are 

15  14  34  32 

temperature  dependent,  such  as  N  /N  and  S  /S  ,  perhaps  increasing 
the  number  of  independent  thermometers  in  bio-organic  data  bases  to 
five  or  six.  This  possibility  should  be  investigated. 

III.  THE  QUESTION  OF  EQUILIBRIUM 

In  the  following  paper  an  attempt  is  made  to  estimate  the  magni¬ 
tudes  of  the  temperature  coefficients  of  the  several  isotope  thermometers 
in  bio-organic  data  bases.  This  can  be  done  using  thermodynamic 
considerations,  if  one  assumes  that  the  organic  material  is  manufactured  in 
eq ui  1  ibriumwithits  i  imnjiir  mmimpH nn. 

but  one  which  does  not  seem  too  unlikely.  For  example,  in  experiments  in 
the  laboratory  the  isotopic  composition  of  CaCC>3  slowly  formed  from 
aqueous  solution,  was  found  to  be  the  same  as  in  the  shells  produced 
by  aquatic  organisms  at  the  same  temperature.  ^  In  any  case,  this 
assumption  can  be  tested  for  plants  and  animals  grown  in  laboratory 


conditions . 


In  the  wild  state,  in  coastal  waters  and  estuaries , 
the  carbon  and  oxygen  isotopic  composition  of  shell  carbonate 
has  been  shown  to  be  in  isotopic  equilibrium  with  bicarbonate 
dissolved  in  the  water. (8)  The  life  span  of  the  shell  fish 
is  about  one  year,  so  it  would  seem  equilibration  in  self¬ 
manufacture  might  take  place  rapidly,  say  in  less  than  or 
about  a  day.  Equilibration  between  C02  gas  and  sea  water 
is  known  to  occur  in  less  than  48  hours, and  between 
l'°2  gas,  water  and  precipitated  CaCO,,  in  less  than  6  days.(10»n) 
In  living  animals,  tritiated  water  is  known  to  equilibrate 
with  blood  serum  of  rats  in  6  days  or  less.(12)  Thus  the 
evidence,  although  meager,  is  that  equilibration  takes  place 
rapidly  compared  with  the  relevant  life  spans. 

It  should  be  pointed  out  that  in  shells  the  carbon  and 
oxygen  are  not  independent  thermometers  because  they  are  tied 
by  the  same  vibrations  in  the  C-0  bond.  But  there  is  also 
organic  material  in  the  shell  matrix  which  could  provide  a 
set  of  multiple  thermometers,  and  furthermore  could  be 
FadiocaTboh  d3LUd.  —  ■ 

The  question  of  the  degree  to  which  equilibrium  obtains 
in  formation  of  cellulose  is  interesting  in  itself.  If 
equilibrium  seems  to  apply,  the  assumption  could  be 
extended  to  computations  of  production  of  other  substances,  if 
equilibrium  does  not  obtain,  and  instead  the  processes  of  isotope 
fractionation  are  kinetic,  then  the  temperature  coefficients  should 
be  even  larger  than  those  estimated  by  assuming  equilibrium. 


(3) 


We  shall  be  computing  partition  functions  at  ^300°K  so  that  the 
harmonic  oscillator  approximation,  is  reasonably  good.  Here  k  is 


the  Boltzman  constant,  t  is  the  absolute  temperature,  v,  is  the  frequency 
til 

of  the  j  appropriate  vibration  of  degeneracy  g^  and  n  is  the  vibrational 
quantum  number,  and  is  the  difference  in  zero  point  energy  between 
the  isotopes.  In  cellulose  the  barriers  preventing  free  internal 
rotation  are  high,  so  that  rotations  deteriorate  into  torsional  os¬ 


cillations  or  bending  vibrations.  In  the  present  treatment  these  will 
be  included  in  as  Par^  of  the  vibrational  partition  function. 

The  stretching  and  bending  vibrations  appropriate  to  cellulose, 
which  is  a  chain  molecule  of  module  (H  -  C  -  Oil),  are  listed  in  Tables 
I  and  II.  Their  frequencies  differ  slightly  for  molecules  containing 
different  isotopes  because  the  isotopic  mass  enters  into  the  reduced 


mass,  y,  of  the  appropriate  oscillator  from  which  the  particular 
frequency  is  computed,  according  to: 


where  K  is  the  force  constant  of  the  particular  vibration.  Substituting  in  (3) 
the  frequencies  from  Table  I,  and  noting  that  the  torsional  frequency, 

920  cm  \  is  3-fold  degenerate,  one  computes  Q(273°K)  =  1.2302  and 
Q(298°K)  =  1.2981  for  the  module  {H  -  C  -  OH}.  In  principle,  the 
D-H  isotope  effect  in  cellulose  involves  two  possible  substitutions:  namely 
substitution  of  a  deuteron  for  the  hydrogen  linked  to  the  carbon  atom  and 


Table  I 


Bond 


Stretching  Vibration 


C  -  H 

i 


i  i 

c  -  c  - 

I  I 


-  0  -  H 


-  C  -  (OH) 

i 


o16  -  o16 


2960  cm 


900  cm 


3680  cm 


(1200  cm'1)* 


1580  cm  ^  B  “  1.44560  cm  ^ 


*  Calculated  from  a  force  constant  taken  as  1/2  of  force  constant 
for  C  »  0,  see  G.  Herzberg,  Infrared  and  Raman  Spectra.  D.  Van 
Nostrand  &  Co.,  Inc.,  New  York,  1945,  Table  51,  pg.  155,  and  see 
also  Table  89. 


Table  II 


Bending  Vibration 

Bending  Frequency 

m  1 

R  -  C  -  OH  • 

920  cm 

R 

?*\ 

R  -  C  -  OH 

920  cm 

R 

H 

1 

R  -  C  -  OH 

920  cm” 1 

H 

R  -  C  -  0  -  H 

| 

(700  cm  )* 

R 

• 

H 

i 

R  -  C  -  OH 

fe  • 

Sr 

(375  cm  *)** 

*  See  Herzberg,  Infrared  and  Raman  Spectra,  pj.  cit.,  page  196,  and 
Tables  51  and  89. 


**  See  Herzberg,  Infrared  and  Raman  Spectra,  0d,  cit..  Tables  118  and 

119. 


secondly  substitution  of  the  hydrogen  attached  to  oxygen.  In  the  case 
of  R  and  R'  attached  to  the  pivotal  carbon,  the  appropriate  reduced  mass 
is  insensitive  to  substitution  of  H  by  D  because  the  masses  of  R  and  R' 
are  almost  infinitely  large  relative  to  H  or  D.  The  partition  functions 
computed  for  the  various  isotopic  substitutions  of  the  cellulose 
module  are  given  in  Table  III. 

Because  -of  the  crudeness  inherent  in  assuming  bond  frequencies 
for  cellulose  to  be  equal  to  those  known  in  rather  small  organic 
molecules,  we  neglect  anharmonic  corrections,  and  complications  such 
as  hydration  in  solution,  effect  of  hydrogen  bond  formation,  etc. 

The  fractionation  ratios  may  be  written  as  follows.  For  the 
oxygen  isotopes,  there  are  three  equilibria  correspond ir.g  to  isotopic 
exchange  between  cellulose  and  each  of  the  three  oxygen-containing 
molecules  involved  in  the  reaction: 


K18(02)  .  <j*(c'i2°18)  Q(0216) 

Q(ch2o16)  q«(016013) 

k18(co  1  .  Q*(CH2 °18)  q(co216  > 

2  Q(ch2o16)  q*(co16o18) 

^  v  .  -  I  .  «ll» 

k18(h20)  -  ^°18)-  q<V16> 

q(ch2o16)  q*(h2o18) 


(5a) 


(5b) 


(5c) 


where  Q*  is  the  partition  function  for  the  heavier  isotope. 


Table  III 


I 

Partition  Functions  for  the  Module  of  Cellulose  (H  -  C  -  OH) 

I 

and  Its  Isotopic  Modifications 


Q(HCOH)  Q*(DCOH)  Q*(HCOD)  Q*(C13H20)  Q*(CH2018) 


1.2302 


40.850 


35.014 


1.4240 


1.5620 


1.2981 


32.695 


28.086 


1.4844 


1.6659 


For  the  hydrogen  isotopes 


two  equilibria  enter  so  that 


Q*(HC0D) 
K(HCOD)  =  _ 

q(ch2o) 


Q(H„0)  Q*(DCOH)  Q(H  0) 

;  K(DCOH)  =  _  _  (6) 

Q*(HD0)  Q(CH20)  Q*(HD0) 


and  for  the  carbon  isotopes,  only  one  equilibrium  enters,  that  with 
C02 ,  so , 


q*(c13h2o) 

q(c12h2o) 


Q(c12o2) 

q*(c13o2) 


(7) 


but  what  the  chemist  actually  measures'  )  is  an  isotopic  fractionation,  a. 
For  example  for  oxygen  the  fractionation  factor  corresponding  to 
equation  5a  is: 


{2[0218]  +  [018016]}  [CH2016] 

{[018016]  +  2[0216]}  tCH2018] 


Considering  the  reaction, 


2t018016] 


(9) 


the  equilibrium  constant  is  4  ,  neglecting  a  tiny  correction  for  isotope  effectf6^ 


4 


fo18o161  2 
to218no216 


] 


(10) 


because  O^8  and  0,,^ »  being  symmetrical,  have  only  half  as  many 

16  18 

rotational  states  as  the  asymmetrical  molecule  0  0  .  Substituting 

(10)  into  (8),  the  fractionation  factor  reduces  to: 


18  Q[o216]1/2  q*[ch2o18] 

3  W2>  qS/5? 


18  16 

Correspondingly  C02  and  C02  ,  being  symmetric,  have  only  half  the 

18  16 

rotational  states  as  C0X  0  ,  so  that 


18....  <!‘co216]1/2<!S!CH2o181 
3  (  2>  V[C0218]l/2  QtCH2016] 


but  for  a(HCOD)  and  a(DCOH)  the  isotope  effect  is  large  so. 


a(HCOD) 


Q*[HCOD]  [2Q[H20]  +  Q[HDO] 

Q[CH20]  [Q[KD0]  +  2Q[D20] 


a(DCOH) 


Q*  [DCOH]  2Q[H20]  +  Q[HD0] 


Q[CH20]  QtHDO]  +  2Q[D20] 


Here,  while  H20  and  D20  are  symmetric,  none  of  the  deuterated  or 

13  12 

hydrogenated  cellulose  modules  are.  Finally  because  C  02  and  C  02 

16  1 R 

have  the  same  number  of  rotational  states,  also  H20  and  H20  ,  the 

fractionation  factors  are 


18  Q(H2016)  Q*(CH  O18) 

a  (H«0)  —  i  r-n-  -I  r 

Q*(H2Oi0)  Q(CH2Oib) 


q(c12o2)  q*(c13h2o) 

q*(c13o2)  q(c12h2o) 


13  12  18  16 

For  the  isotope  ratios  D/H,  C  fZ  ,  and  0  /0  ,  we  may  now  compute 


the  partition  function  ratios  CQ*/Q) 


^  q.  The  partition  function  ratios 


for  (Q*/Q)oxygen  are  taken  from  Urey  (Ref.  6),  those  for  (H2°^liq  are 

(13) 

taken  from  computations  of  Bottingav  and  those  for  (C02)gag  from 
Bottinga^3^  and  for  HDO  and  D20  from  Urey^ 

The  ratios  appropriate  for  each  isotope  exchange  are  listed  in 
Tables  IV,  V  and  VI.  By  substituting  these  in  equations  (11)  through 
(15)  the  fractionation  factors  have  been  calculated  for  273°K  and  298°K. 
The  corresponding  temperature  coefficients  are  listed  in  Table  VII. 

•  1  o  in  _ 

The  temperature  coefficient  so  calculated  for  C  /C  of  0.36  ppt/  C 

13  12 

agrees  with  the  coefficient  measured  for  C  /C  in  plant  and  animal 

(3  4) 

material  by  Sackett  and  his  co-workers,  *  so  that  the  assumption  of 
equilibrium  appears  to  be  more  or  less  valid  and  one  can  hope  that  the 
like  calculations  for  oxygen  and  hydrogen  may  be  meaningful. 


From  Ref.  6,  at  equilibrium  the  temperature  coemgtent—frf  - 
HD0/H20  is -312%o  per  degree  C,  namely  much  larger  than  that  of  the 
(D/H)  ratio  calculated  for  cellulose,  so  that  measuring  it  in  a 
geographical  distribution  of  the  same  kind  of  trees  would  measure  the 


— ±fc>- 


Table  IV 


Partition  Function 

273°K 

298°K 

Comment 

Q*[DC0II]/Q[HC0H] 

33.206 

25.187 

Table  III 

Q*  [HCOD] / Q [ HCO  H] 

28.462 

21.636 

<Q*[D2OJ8as«[H2°]gils)1/2 

16.503 

12.543 

H.  C.  Urey  (Ref.  6)* 

Q*[HD0]g3s/Q[H20]gas 

32.7400 

24.9460 

*  The  partition  function  ratio  for  liquid  is  obtained  by  multiplying  the 
partition  function  ratio  for  gas  by  the  D/H  ratio  in  liquid  divided  by 
that  for  gas.  It  is  taken  here  as  1.104  at  273°K,  and  as  1.074  at 
298°K,  where  the  numbers  in  parentheses  refer  to  partition  of  tritium, 
see  Ref.  22. 


Table  V 


Partition  Function 

273°K 

298°K 

Comment 

q*[c13h2o]/q[c12h2o] 

1.1575 

1.1435 

Table  III 

In  Q*[C1302]/Q[C1202]  (gas) 

0.19732 

0.17558 

9 

Y.  Bottinga 

Q*[C1302]/Q[C1202]  (gas) 

1.2181 

1.1919 

Table  VI 


Partition  Function 

273°K 

298°K 

Source 

• 

q*[ch2o18]/q[ch2o!16 

1.2697 

1.2833 

This  paper 

[Q**(028)/q(026)]g/f 

1.0923 

1.0818 

Urey  (Ref.  6) 

1/2  ln{Q*[C028]gas/Q[C026]gas} 

0.12530 

0.11108 

Bottinga,  (Ref. 

13) 

(Q*[C028]/Q[C026]}1/2gas 

1.1336 

1.1175 

ln{Q*[H2018]gas/Q[H20l6]gas} 

Q*[H2018]/Q[H2016]gas 

0.06822 

1.0706 

0.06164 

1.0635 

Bottinga,  (Ref.  13) 
See  comment  below 

Comment:  The  partition  function  ratio  for  liquid  water  is  obtained  by 
multiplying  the  partition  function  ratio  for  water  vapor  with 
the  0l8/0lb  ratio  in  liquid  water  divided  by  the  ratio  in  the 
vapor.  This  correction  is  taken  as  1.01150  at  273°  and  as 
1.00930  at  298  K;  (Ref.  13,  page  806) 


Tabic  VII 


\ 


• 

273°K 

298°K 

Temperature 

Coefficient 

a18(°2) 

1.162A 

1.1863 

+0.96 

o ,  o„ 

/  per  C 
oo 

a18(C02) 

1.1200 

1.148A 

+1.1A 

o  ,  o_ 

Lr,  Per  C 

a18(H20) 

1.1725 

1.1956 

+0  92 

o ,  o_ 

/  per  C 
oo  v 

13 

a 

0.9503 

0.959A 

+0.36 

o ,  o_ 

/  per  C 

oo  r 

a  (HCOD) 

1.5881 

1.5983 

+0.  A 

o ,  on 

/  per  C 
oo  * 

a  (DCOH) 

1.8095 

1.8606 

+2.0 

o .  o„ 

/  per  C 
oo  * 

The  computed  temperature  coefficients  listed  above  agree 
very  well  with  the  \alue  0.35  ppt  per  C  measured  for  Cl3/ 
Cl2  by  Degens  et.  al.,  Ref.  4,  and  with  the  value  3.0  ppt 
per  °C  measured  for  D/H  by  W.  E.  Schiegl,  Science  175,  512- 
513,  (1972),  "P^uteriura  Content  of  Peat  as  a  Paleo-Cliraate 
Recorder. " 


distribution  in  D/H  ratios  of  rain  water,  and  measuring  it  in  a  chronological 
sequence  would  measure  time  dependence  of  the  rain  water. 

Our  naive  assumptions  are  mainly  as  follows.  We  have  assumed 
that  cellulose  is  formed  in  solution,  that  in  cellulose  rotations  have 
degenerated  into  torsional  vibration,  that  its  vibrations  are  like  those 
in  small  organic  molecules,  and  that  manufacture  of  cellulose  in  plants 
is  an  almost  equilibrium  process!'^^^ 

The  intent  of  the  present  computation  is  to  show  that  all  three 
of  the  isotope  ratios  D/H,  C13/C12,  and  01S/016  nay  well  be  thermometers, 
and  that  their  temperature  coefficients  may  have  measurable  magnitudes 
in  bio-organic  plant  residues,  and  so  may  give  information  about 
climate  variation  in  the  past. 

For  example,  we  may  measure  these  ratios  in  a  chronological 
sequence  of  tree  rings.  The  results  of  Parker (l5>  are  that  (1)  different 
plant  species  growing  side  by  side  in  a  marine  estuary  can  have  different 
percentages  of  C13  in  total  carbon  content,  and  (2)  individuals  of  the 
same  species  have  about  the  same  percentage  of  C13.  So  in  evaluating 
past  climates  from  isotope  thermometers  we  should  measure  changes  in 
a  single  species  of  tree  and  we  should  measure  a  chronological  sequence 
consisting  of  several  individuals,  of  the  same  kind  of  tree,  having  over- 
lapping  life  spans.  By  measuring  trees  of  overlapping  life  spans,  and 
trees  that  grow  In  sparse  environments  such  as  mountain  slopes,  one  may 
avoid  spurious  effects  caused  by  the  relatively  greater  removal  of 
C02  from  air  near  the  ground  by  abundant  low  growing  plants. 


From  Ref.  6,  the  temperature  coefficient  of  oxygen  in  rain  water 
could  be  as  large  as  -0.36%o  per  degree  C.  Also  important  is  the,  fact  that 
CO2  mixes  rapidly  through  the  entire  global  atmosphere  in  a  time  of  ^5 
years,  and  with  the  oceans^1^  ^  in  about  15  years,  so  that  at  any  time 
the  atmospheric  carbon  isotope  ratios  reflect  the  temperature-dependent 
isotope  separation  in  the  sea,  which  has  a  coefficient  of  -0.21  ppt/°C 
for  oxygen  (Ref.  10),  and  -0.11  ppt/°C  for  carbon  (Ref.  11). 

One  may  ask  whether  there  is  isotope  exchange  between  the  hydrogens 

bound  in  old  cellulose  and  in  new  sap.  As  shown  in  Table  I,  both 

hydrogens  are  bound  very  tightly,  even  more  so  than  the  C-C  bonds,  so  that, 

knowing  that  radio-carbon  dates  in  tree  cellulose  mainly  agree  with  the 

tree  ring  dates,  meaning  that  in  heart  wood  carbon  does  not  exchange,  we 

can  hope  the  sane  to  be  true  of  hydrogen  in  heart  wood.  But  this  point 

should  be  tested.  When  cellulose  is  prepared  in  a  thin  layer  or  finely 

(21) 

ground,  and  repeatedly  wetted  and  dried,  there  is  some  exchange, 
but  these  conditions  are  quite  different  from  cellulose  in  sound  heart 
wood . 
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PRINCIPLE  OF  THE  METHOD 

Preliminary  measurements  indicate  the  feasibility  of  using  fast- 

neutron  reactors  to  measure  0  /0  isotope  ratios  with  an  accuracy  of 

4 

1/10  ,  by  neutron  activation.  A  fast  reactor  is  necessary  because  its 

13  2 

total  neutron  flux  is  10  /cm  sec  or  more  and  because  its  ultra  fast 

1  ft 

neutron  flux,  F  ,  and  slow  neutron  flux,  F  ,  are  comparable.  The  0  , 

u  s 

with  a  0.2/  natural  abundance,  is  activated  by  thermal  neutrons  via  the 

3 

16 

(n,v)  reaction,  while  the  0  is  activated  via  the  (n,p)  reaction  by 

I 

M 

ultra-high  energy  neutrons  (>10  MeV)  in  the  high  energy  tail  of  the 

I 

fission  neutron  spectrum. 

Both  oxygen  activities  are  measured  simultaneously  follow¬ 
ing  activation,  and  ideally  should  have  approximately  equal  numbers  of 
decays  or  counts,  for  optimum  determination  of  the  ratio  R  =  018/016. 

This  equality  can  be  achieved  because  the  thermally  activated  isotope 

18 

happens  to  be  the  rare  one  (0  ),  so  that  its  low  abundance  is  compen¬ 

sated  by  a  flux  of  thermal  neutrons  larger  than  the  fast  flux. 

The  characteristics  of  these  activities  are  as  follows: 

(a)  0^8(n,p)N18;  threshold  neutron  energy  «=  10.25  MeV, 

°l6(n»P)  e  4°  mb,  t1g  (mean  life)  =  10.3  sec 
■=  6.14  MeV  (68%) 

(b)  0  (n,Y)0^;  formed  by  thermal  neutron  capture, 

°18^Y,n^  =  0*22  mb»  fci8  (raean  life)  *  41.8  sec 
Eyl  ■  1.37  MeV  (60%),  E  2  =  200  KeV  (100%) 


For  long  irradiation  times. >,  60  sec,  equilibrium  is  achieved,  and  the 


ratio  RprtOf  the  total  decays!  detected  is  given  by 

/ 


Fs  °18  fc18  e18  = 

F  a.,  t..,  e., 

u  16  16  16 


—3  18 

where  a  =  2  x  10  (abundance  of  0  );  e^g  and  are  efficiencies  of 

detection  of  the  photo  peak  for  the  respective  gamma  rays,  with  a  ratio 


e18^e16  ^or  a  ^  x  ^  Nal  detector  crystal  (Harshaw  Chemical  Company 

catalog,  "Scintillation  Phosphors,"  1962).  With  these  figures 


*EQ  ’  4'6  •  10'‘  VF„ 

For  short  irradiations  $  7  sec,  the  ratio  R  of  accumulated  counts  is 

9  s 

Rs  c  t16  ^0^8  =  1,1  ’  10  REQ 

The  optimal  ratios  of  ultra  fast  to  slow  flux  for  and  Rg  unity  are 

3  4 

respectively  2.2  •  10  and  10  ,  These  flux  ratios  are  approximately 
obtainable  in  the  TRIGA  fast  reactor  at  Gulf  General  Atomic  Corporation 
(GGA)  in  La  Jolla,  California.  Accordingly,  one  of  us,  Dr.  Richard 
Lukens,  manager  of  activation  analysis  of  GGA,  verified  the  simultaneous 
and  observable  production  of  both  oxygen  activities  by  neutron  activation. 
These  results,  shown  in  Figure  1,  display  clearly  both  the  200  KeV  and 
the  6  MeV  gamma  rays  in  both  wood  and  water.  The  conditions  of  his 
experiment  allow  a  great  deal  of  improvement  to  be  made  for  the  accurate 
measurement  of  R. 

Based  on  the  successful  proof,  achieved  in  the  GAA  experiment, 
that  R  can  be  measured  by  neutron  activation,  in  the  present  report 
we  analyze  steps  for  improvement  of  measurement.  The  accuracy  which  is 


COUNTS 


4 

ultimately  desired  is  1/10  ,  this  being  the  accuracy  advertised  for 
mass  spectrometers  (MS),  however,  the  MS  measurements  are  made  on  C^, 
prepared  and  purified  by  several  successive  steps  in  some  of  which 
isotope  fractionation  is  an  uncomfortable  possibility. 

If  the  requisite  accuracy  can  be  achieved,  measuring  R  by  neutron 
activation  may  offer  the  advantages  of  checking  on  the  mass  spectrograph 
results  by  an- alternate  method  and  avoiding  possible  chemical  and  mechani¬ 
cal  fractionation.  Non-destructive  measurement  of  the  sample  may  also 
be  .a  possibility.  So  this  new  method  might  be  very  valuable.  Once  de¬ 
veloped  and  standardized,  a  number  of  practical  applications  are  expected 
to  follow. 

METHOD  AND  NUMERICAL  ESTIMATES 

18  16 

In  designing  a  method  of  measuring  the  0  /0X  ratio  by  neutron 

activation,  there  are  three  main  problems.  The  first  is  to  find  a  distinct 

means  of  identification  for  each  isotope,  then  come  the  two  related 

problems  of  interfering  activities  and  of  statistical  accuracy.  The 

16  18 

identification  tags  for  the  0  and  0  isotopes  could  consist  of  three 

19 

different  gamma  rays  produced  by  neutron  activation,  two  for  0  (from 
the  018(n,Y)  reaction)  and  one  for  N16  (from  the  016(n,p)  reaction). 

IQ 

Figure  2  gives  the  detailed  decay  schemes  for  0  and  N  [3],  Here  we 
concentrate  on  the  6.143  MeV(N^"8)  and  0.200  MeV(O^)  y-rays.  The  reason 
for  this  is  that  the  1.37  MeV  gamma  ray  is  completely  obscured  in  water 
by  traces  of  sodium  and  in  wood  by  sodium  and  other  elements.  (See 
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Figure  2  Decay  schemes  for  190  and 
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Figure  1)  By  chance,  neutron  activation  of  sodium,  which  is  a  persistent 
impurity,  produces  a  gamma  ray  at  1.36  MeV,  unresolvable  from  the  O18 
radiation  at  1.37  MeV.  This  obscuration  was  noted  in  an  initial  experi¬ 
ment  made  with  the  medical  reactor  at  Brookhaven  National  Laboratory  (BNL) , 
with  Dr.  Hobart  Kraner,  (described  in  RDA-MR-031-NMO,  March  31,  1972). 

As  already  noted,  measuring  R  requires  both  a  copious  ultra-fast  neutron 

flux  (E  >  10,2  MeV)  and  a  thermal  flux  whose  total  should  be  ~1013/cm2  sec 

or  more,  so  that  the  required  statistical  accuracy  can  be  achieved  in  a 

reasonable  time.  The  point  we  make  now  is  that  the  identification  and 
16  18 

measurement  of  0  and  0  by  their  respective  y  peaks  using  Nal  crystal 
detectors  has  been  shown  to  be  possible  in  the  GGA  experiment. 

We  shall  concentrate  on  wood  because  of  its  importance  to  the  study 
of  climatology;  namely  the  variation  of  isotopic  ratios  of  such  elements 

as  C,  0,  H,  N,  S,  - ,  in  tree  rings,  is  expected  to  depend  on  the  ambient 

air  temperature  when  the  tree  ring  was  formed,  and  therefore  to  have  re¬ 
corded  variations  of  climate. 

The  major  elements  of  interest  in  wood  are,  of  course,  C,  II,  and  0. 

We  consider  interference  in  measurement  of  R  caused  by  activation  of  trace 
elements  in  wood,  such  as  F,  Mn,  N,  S,  K,  which  may  vary  in  abundance 
from  one  species  of  tree  to  another.  However,  interference  caused  by 
trace  elements  can  be  minimized  (except  for  sodium).  For  example,  con¬ 
sider  Mn55  activated  by  thermal  neutrons  with  properties  as  follows: 

Mn55(n,y)Mn56,  thermal  neutron  capture;  o(n,y)  =  13.3b; 

^1/2  18  ^*6h;  =  0.845  MeV(100%);  abundance  of  manganese 

in  wood  a  1.5  x  10~*. 


Since  0  must  be  determined  by  measuring  the  200  KeV  radiation,  there 
will  be  significant  interference  at  ~200  KeV  caused  by  Compton  scattering 
of  the  0.845  MeV  y-ray  into  the  200  KeV  region.  Indeed,  in  the  prelimi¬ 
nary  experiment  at  BNL,  the  200  KeV  gamma  peak  in  irradiated  wood  was 
completely  masked  by  the  high  manganese  content.  But  in  the  BNL  experiment 
the  irradiation  time  was  long  (90  sec)  trad  there  was  considerable  delay 
flfter  extraction  of  the  wood  sample  from  the  reactor  before  counting 

the  200  KeV  gamma,  both  of  which  factors  enhanced  the  Mn  activity  rela- 
18 

tive  to  0  .  A  third  factor  was  that  the  Nal  detector  was  3"  x  3"  so 

that  the  signal  to  noise  ratio  was  poor  at  200  KeV. 

Two  of  these  parameters  were  improved  in  the  syccessful  GGA  experi¬ 
ment.  Namely  the  irradiation  time  was  decreased  to  30  sec  (still  not 
the  optimal  time  of  ~10  sec) ,  and  the  delay  before  counting  was  almost 
minimized.  However,  the  Nal  detector  still  had  a  too  great  thickness 
of  3  in.  The  efficiency  for  detection  of  the  Mn  gamma  ray  of  845  KeV 
which  has  been  Compton-scattered  into  the  200  KeV  region  can  be  de¬ 
creased  by  about  a  factor  of  10  in  a  thin  Nal  crystal,  without  de¬ 
creasing  the  efficiency  of  measurement  of  the  200  KeV  y  ray  emitted  by 
oxygen.  A  thin  Nal  crystal  will  be  used  in  the  next  experiment  at  GGA 
to  achieve  this  enhancement. 

The  GGA  reactor  can  be  pulsed  at  ~10^2  neuts/cra2  sec  flux  for 
3 

x  10  sec  for  a  total  flux  of  7  •  10^  neuts/cm2,  as  compared  with 
the  total  flux  of  ~1013  neuts/cm2  used  in  obtaining  the  data  of  Figure 
l.l  i  Two  advantages  may  result  from  irradiation  by  pulsing,  namely, 
the  short  irradiation  time  minimizes  the  effect  of  impurities  which  have 
long  lived  activities,  and  secondly,  the  fast  neutron  flux  is  enhanced, 


which  increases  the  016(n,p)N36  activity  relative  to  the  O1^  activity. 

Pulsed  activation  will  be  explored  as  a  method  to  improve  measurement  of  R. 

Interfering  Activations:  Radio-active  can  be  produced  by  other 

16 

reactions  than  activation  of  0  ,  causing  the  measured  value  of  R  to  be 

too  low.  These  reactions  are  listed  in  Table  I. 

Table  I 

ACTIVATION  REACTIONS  PRODUCING  RADIOACTIVE  N16 


Production  Cross  section 

Isotopic 

Abundance 

Threshold 

Energy 

w19,  ...16 

(n,a)N 

2  x  10-27  cm2  at  5-6’ MoV 

100% 

1.5  MeV 

JL7,  ,...16 
gO  (n,d)N 

- 

.037% 

11.5  MeV 

-18,  ....16 
gO  (n,t)N 

- 

.204% 

13.3  MeV 

..15,  .16 

?N  (n,y)N 

-29  2 

~2  x  10  cm 

.37% 

0 

We  must  ask  whether  the  radio-N^  produced  from  any  ^*f  these  reactions 

16  -4 

by  irradiation  of  wood  in  a  reactor  can  increase  N  by  more  than  10 
oi  the  production  from  0^(n,y)N^. 

According  to  Figure  1,  the  ratio  R  measured  in  the  GGA  experiment  is 
about  10  (instead  of  the  desired  value  of  '“l),  but  if  the  irradiation 
time  which  was  30  sec  had  been  shortened  to  the  desired  time  of  10  sec, 


and  if  the  waiting  time  before  counting  had  been  decreased  from  30  to 
10  sec,  the  ratio  R  would  have  been  close  to  unity,  so  that  F  /F  ~  10^ 
in  the  TRIGA  reactor  for  fast  neutrons  of  $  10  MeV.  We  may  guess  the 
relative  flux  of  neutrons  at  the  intermediate  energy  of  5  MeV  to  be  30 
times  greater.  Then  the  N16  activity  produced  from  F19(n,a)N16  relative 
to  that  from  016(n,p)N16  is  approximately  less  than  or  equal  to 


where  the  first  bracket  is  the  ratio  of  abundances  of  fluorine  and  oxygen, 
the  second  is  the  ratio  of  cross  sections 


Thus,  the  fluorine  content  of  wood  should 
not  be  a  nuisance.  For  experimental  evidence  of  the  fluorine  con¬ 

tent  we  may  look  for  the  reaction  F19(n,Y)F20,  T±/2  =  11  sec,  Ey  =  1.63  MeV, 
which  should  appear  at  channel  ~52  in  Figure  1  as  a  peak  in  wood  (but  not 
in  water),  with  an  intensity  of  10  relative  to  N16,  however,  no  peak  is 
there.  In  fact,  wood  and  water  are  not  perceptibly  different.  The 

answer  may  be  that  fluorine  is  less  than  10  ppm  in  this  particular  kind 
of  wood. 


The  basic  formula  for  cellulose  is  (cgKjo°5^x  Present  in  wood  to  (Fig.  3) 
about  70%.  Lignin,  with  composition  (C^H^O^,  is  present  to  about 
25%,  and  resins,  (C^HgO^,  to  a  few  percent.  The  elemental  impurities 
of  interest  are  present  in  amounts  typically  as  follows:  N  at  ~0.2%, 

10  ppm  for  S,  100  ppm  for  K,  10  ppm  for  F,  20  ppm  for  AX,  40  ppm  for  Na, 
and  150  ppm  for  Mn. 

Because  of  the  high  energy  threshold  for  g017(n,d)N16,  the  ratio  of 
M16  ,  , 

N  it  produces  is  given  by 
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Figure  3 


g(mb)  3.7  x  10 
40  0.5 


1  =  2  x  10 


o(mb) 


so  that  for  a  reasonable  cross-section  of  a  few  mb  or  less,  this  reaction 

18  16 

is  negligible.  The  ratio  for  gO  (n,t)N  ,  with  a  similarly  high  energy 
threshold,  is  given  by 


c(mb)  (2  x  10~3) 
40  |  0.5  i 


{1} 


o(nb) 


and  even  with  a  grossly  overestimated  cross-section  of  e*l  mb  the  contri¬ 
bution  of  this  reaction  would  be  tolerable. 

Next,  let  us  consider  ^N^ (n,y)N^,  produced  by  thermal  neutrons. 
Nitrogen  occurs  in  the  nucleic  acid  attached  to  lignin,  to  an  amount 
which  could  be  as  large  as  0.2%  of  the  mass  of  wood.  The  activity 
it  produces  is  given  by,  relative  to  the  activity  from  O1^ 

l-x  10 — -)  104  S  7  •  10~5 

4  x  10"26j 


(3.7  x 


-3 

10 


(2 


x  10~3)] 


0.5 


where  the  first  bracket  is  the  isotopic  abundance  of  N15,  the  second  is 
the  ratio  of  elemental  abundance  of  N  and  0  in  wood,  the  third  is  the 
cross  section  ratio  and  the  fourth  is  the  appropriate  ratio  of  neutron 
fluxes. 

The  only  contributor  to  the  200  keV  019  peak  is  F19  from  the 
19  19 

F  (n,p)0  reaction.  The  ratio  of  its  contribution  to  the 
0^®(n,y)0^9  contribution  is  given  by 


oo 


/_ . . ?■  \  /10  %  /  30  mb  >  /10  N  "*< 

v  To-;  ( - )  ( - )  ( - ) 

2.2  x  10  0.5  0.22  mb  1 


-3 

l  x  10  , 


making  it  a  significant  contributor  to  the  error.  If  however  one 
utilizes  the  1.63  MeV  gamma  ray  from  the  F19(n,y)F20,  tt  =11  sec 
reaction,  one  can  relate  the  number  of  counts  in  the  1.63  MeV  peak  to 
the  contribution  of  F  to  the  200  keV  peak  and  subtract  it.  The  ratio 
of  the  F  (n,p)0  reaction  to  the  F*9(n,y)F2^  reaction  need  only  be  known 
to  1/10  (from  an  independent  measurement  in  the  system)  to  reduce  the 
F  (n,p)0  contribution  ratio  to  10  .  It  is  possible  to  measure  the 
fluorine  reaction  ratio  to  at  least  1/10  thereby  making  the  fluorine 
(n,p)  contribution  insignificant. 

So  it  seems  that  measurement  of  0  /0^  in  wood  may  possibly  be  made 

non-destructively  by  neutron  activation. 

The  statistical  error  in  the  measurement  of  the  0^/0^  ratio  is 
the  error  in  the  total  number  of  counts  under  each  of  the  gamma  peaks 
of  interest.  This  error  for  107  counts  is  ~l/2000. 

The  problem  that  must  be  solved,  in  th?.  nuclear  activation  method  for 
determining  the  0  /O1^  ratio  with  an  accuracy  of  «3/10^,  is  recording 

high  count  rates  without  encountering  unacceptable  losses  in  count  rate. 
This  problem  arises  because  N16  has  a  half-life  of  7.2  sec  so  that  the 
effective  data  collection  time  is  »15  sec.  To  collect  107  counts  in  this 
time  requires  an  initial  count  rate  of «2-3  x  106  sec  in  the  peak  of 
interest.  The  technology  of  counting  at  megacycle  rates  exists,  namely 
using  fast,  solid  state  discriminators,  and  this  problem  we  are  next 
attacking. 
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Determination  of  Isotopic  Abundances  of  C,  N,  and  0  in  Organic 
Samples  Using  Nuclear  Techniques  * 


K.  W.  Jones  and  H.  W.  Kraner 
Brookhaven  National  Laboratory,  Upton,  New  York  11973 


ABSTRACT 

The  determination  of  isotopic  abundances  of  C,  N, 

and  0  in  organic  samples  by  techniques  involving  the  use 

of  charged  particle  beams  from  a  Van  de  Graaff  accelerator 

is  discussed.  The  accuracy  and  convenience  of  the  method 

are  considered  and  an  example  of  a  method  for  measuring 
16.  18 

0-  0  ratios  is  given. 


*Work  performed  tinder  the  auspices  of  the  U.S.  Atomic  Energy  Commission. 


I.  Introduction 


Radioactive  isotopes  have  been  used  for  many  years  as  tracers  in 
biological  and  medical  experiments.  The  activity,  of  course,  is  used  as 
a  tag  to  trace  the  progress  of  the  isotope  through  the  sample.  While 
this  technique  is  very  powerful  it  does  have  the  disadvantage  in  some 
cases  that  the  dosage  can  produce  undesirable  radiation  damage  effects 
in  the  sample  or  patient  and  that  great  care  in  handling  of  the  radioactive 
isotope  is  required. 

A  possible  alternative  approach  would  be  to  use  compounds 

which  are  enriched  in  some  isotope  present  in  the  normal  element  in 

only  a  small  percentage . This  might  be  particularly  advantageous  for 

carbon,  oxygen,  and  nitrogen  where  the  radioactive  isotopes  either  have 

14 

relatively  short  lifetimes,  or  in  the  case  of  C,  a  long  lifetime  with 
a  soft  fl-ray,  and  thus  requ ires  relatively  large  amounts  of  activity  and 
special  care  in  the  processing  and  counting  techniques.  One  of  the  major 
and  limiting  difficulties  in  the  use  of  the  stable  isotope  technique  is 
to  achieve  adequate  sensitivity  for  the  detection  of  small  changes  in 
the  isotopic  abundance  of  each  element. 

Many  approaches  to  this  problem  can  be  envisaged,  but  here  the 
discussion  will  be  limited  to  the  determination  of  C,  N,  and  0  isotopic 
abundance  ratios  in  organic  samples  by  nuclear  techniques  which 
involve  the  use  of  charged  particles  easily  available  from  a  4.0  MV 
Van  de  Graaff  accelerator.  This  type  of  accelerator  is  chosen  for  con¬ 
sideration  since  many  are  available  for  use  and  the  operating  costs  are 

modest.  Naturally,  accelerators  capable  of  producing  higher  energy 

» 

particles  could  also  be  used.  The  following  sections  consider  the 


suitability  and  ease  of  application  of  the  various  nuclear  reactions 
which  can  be  induced  by  these  particle  beams.  Generally,  the  experimental 
techniques  are  fairly  standard  in  nuclear  physics,  but  it  should  be 
emphasized  that  some  care  and  skill  will  be  necessary  for  their  successful 
application,  i.e.  the  user  is  not  completely  free  from  responsible 
knowledge  of  nuclear  techniques. 

II .  General  Considerations 

The  concepts  from  nuclear  physics  which  are  needed  for  the 

determination  of  isotopic  abundances  are  rather  simple  and  have  been 

discussed  many  times  in  the  literature.  For  cbnvenience  and  clarity  a 

brief  discussion  will  be  given  here. 

The  identification  of  a  particular  element  is  made  by  bombarding 

a  target  with  a  beam  of  particles  and  then  measuring  the  energy  of  the 

particle  scattered  by  the  target  or  the  energy  of  some  reaction  product. 

The  kinematics  of  nuclear  reactions  and  scattering  are  well  discussed 

(2) 

by  Marion  and  Young. 

The  energy  release  in  a  nuclear  reaction  is  denoted 

by  Q  and  is  defined  by: 

Q  =  (M:  +  M2  -  M3  -  M4)c2  1) 

where  M^,  M2>  M3>  are  the  masses  of  the  incident;  target.  Hght 
and  heavy  reaction  products,  respectively.  If  and  Q=0  we 

have  elastic  scattering,  if  Q<0  inelastic  scattering,  and  for  a 

reaction. 

The  energy  of  the  light  product  particle  is  the  quantity  which  is 

e 

measured  and  for  elastic  scattering  it  is  given  by: 


2) 


Ei  ^  '  (cos;  +  [(M^M^2-  sin2;]172)2 

Ec  =  (Mt+H2)2 

where  E  is  the  incident  energy,  E.  the  final  energy,  and  ;  the  angle 
o  A 

of  observation  relative  to  the  incident  beam.  Note  that  the  value  of 

E  /E  called  the  scattering  ratio,  is  independent  of  the  incident 

1  ° 

energy.  The  energy  of  the  light  particle  from  a  reaction  is  given  by: 

/  21 12  2 

Ej^E^M2  {(MjM^) 1/2  cos;  +  (MM40/Eq  +  sin  ;)  }  3) 

where  M-j+M^  and  the  energies  are  defined  above.  The  measurement 

of  the  final  energy  of  the  light  particle  at  an  angle  ;  thus  allows  the 
determination  of  the  mass  of  the  heavy  target  particle,  and,  in  particular, 
for  the  case  of  isotopes  of  C,  N,  and  0  it  uniquely  identifies  the  element 
and  isotope. 

The  number  of  atoms  of  a  particular  target  species  present  is  pro¬ 
portional  to  the  number  of  scattered  particles  and  is  found  by  the  following 
relation: 

T  Of  N(Ej ,  ;)/ (Ng  •  c<G)  •  dfi)  *) 

2 

where  T  is  the  number  of  target  atoms  per  cm  , .  N(Ejt;)'  is  the  number  of 

detected  particles  with  energy  E^  at  an  angle  ;,  N^  is  the  number  of  incident 

2 

particles,  c(6)  is  the  differential  cross  section  in  cm  /sr,  that  is,  the 
probability  for  the  reaction  or  scattering  of  interest,  and  dQ  is  the  solid 
angle  subtended  by  the  detector  at' the  target.  The  absolute  accuracy  of  the 
determination  of  T  is  mainly  limited  by  our  knowledge  of  c(6)  which  at  energies 
and  angles  of  interest  for  C,  N,  and  0  measurements  will  be  typically  a 
few  per  cent.  This  means  that  absolute  determinations  of  abundances 
will  be  limited  to  that  accuracy.  However,  if  we  consider  the  measurement 
of  changes  in  the  abundance  of  an  isotope  relative  to  the  composition  of 
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a  standard  sample,  the  accuracy  can  be  much  better;  with  a  practical 
limit  of  ~  0. 17.. 

Abundance  ratios  can  be  found  most  accurately  by  using  a  single 
detector  to  measure  the  energies  of  particles  emitted  at  a  given  angle. 

In  that  case  the  isotope  is  identified  by  the  energy  and  type  of  particle 
and  the  abundance  ratio  is  just: 

$)]•  [ a2  ( e)<7^  (  ]  5) 

If  two  detectors  are  used,  for  example  to  measure  two  different  types 
of  particles,  then  a  ratio  of  solid  angles  must  also  be  included  as 
well  as  possible  differences  in  counter  ef  f  ici'encies .  If  a  standard 
sample  is  used  for  calibration  T^/T 2  Is  known  and  R  can  be  used  to 
determine  the  cross  section  ratio. 

c2  (  6)  /  cx  (  S)=R  •  N2  (Ej  ,  f)  /Nl  (Ej  ,  I)  .  6) 

The  ratio  of  the  absolute  abundances  of  the  two  isotopes  in  the 
unknown  is  then  found  by  substituting  the  ratio  of  cross  sections 
given  in  equation  6  into  equation  5.  The  fractional  abundance 
of  each  isotope  is  then  found  by  remembering  that  the  sum  of  the 
fractional  abundances  is  equal  to  unity. 

Conceptually,  to  this  point,  the  determination  is  trival.  In 
practice  care  in  choosing  detectors  and  obtaining  good  energy  resolution 
will  be  required.  The  most  demanding  task  is  the  use  of  target  techniques 
to  make  thin  targets  able  to  withstand  high  beam  currents  and  the  use 
of  good  vacuum  techniques  to  prevent  contaminants  from  pump  oil  or  other 
vapors  in  the  chambers  from  depositing  on  the  target  and  introducing 
errors.  Such  techniques  are  well  known  to  the  nuclear  physicist;  however 


the  degree  of  ultimate  application  may  well  hinge  on  the  ease  and 
reliability  with  which  targets  can  be  obtained.  One  suggested  method 
to  obtain  oxygen  targets  will  be  discussed  in  a  following  section. 

One  experimental  point  that  should  also  be  mentioned  specifically 
is  the  question  of  target  heating.  Most  organic  materials  cannot 
withstand  the  heating  effects  of  several  microampere  currents  of 
heavy  charged  particles.  Special  care  must  be  given  to  the  prepara¬ 
tion  of  thin  targets  and  to  the  use  of  techniques,  such  as  target 
wobblers,  which  permit  the  use  of  reasonable  beam  currents.  The 
details  will  depend  on  the  compound  of  interest.  We  will  not  dis¬ 
cuss  here  types  of  particle  detectors  appropriate  for  these  experi¬ 
ments  or  details  of  experimental  hardware,  but  refer  the  reader  to 
the  paper  of  Amsel  et  al.  which  gives  an  introductory  survey. 

The  application  of  these  nuclear  techniques  to  qualitative 
analysis  is  by  no  means  new,  but  the  main  emphasis  has  been  on  its 
use  for  trace  element . identification.  A  rather  general  discussion 
of  the  applications  of  elastic  scattering  techniques  has  been  given 
by  Cohen  and  Moyer ^  and  Jolly,  Gruhn,  and  Maggiore.^  A  detailed 
discussion  of  the  use  of  nuclear  reactions  in  measurement  of  light 
element  concentrations  is  given  by  Amsel  et  al.^ 

The  interaction  of  the  standard  beams  of  1H,  2H,  3H,  3He,  and 

4 

He  with  C,  N,  and  0  isotopes  in  the  energy  region  below  4  MeV  have 
been  studied  intensively  for  many  years  and  a  large  literature  exists 
on  the  O-value  and  cross  sections  for  the  reactions.  The  Q-values  are 
tabulated  in  Table  I  and  were  taken  from  the  KAS-NRC  1960  Nuclear  Data 
Tables  and  1971  mass  tables  of  Wapstra  and  Gove . ^  While  later  and 


and  more  accurate  mass  tables  exist,  these  Q-values  are  given  in  a 
convenient  form  in  the  1960  Nuclear  Data  Tables  and  are  more  than 
adequate  for  the  present  discussion. 

Summaries  of  experimental  cross  section  results  have  been 
given  by  Jarmie  and  Seagrave^  and  Kim,  Milner,  and  McGowan.  ^ 

A  very  general  survey  of  all  experimental  results  is  given  by 
Ajzenberg-Selove  and  Lauritsen^  on  A  =  11  and  12  nuclei  and  by 
Ajzenberg-Selove^10^  for  A  =  13-15,  and  A  =  16,17. 

The  rare  isotopes  of  C,  N,  and  0  are  13C,  15N,  1?0,  and  180 
with  normal  abundances  of  1.11,  0.37,  0.037,  and  0.2047.,  respectively. 
Presumably  in  most  cases  the  uptake  of  the  enriched  tracer  compound 
will  not  change  the  isotope  ratios  by  a  very  large  amount  and  it 
will  be  necessary  to  measure  the  changes  with  modest  precision. 

The  precision  of  the  determination  will  be  limited  by  the  counting 
statistics  and  ultimately  the  cross  .section  for  the  process  that 
identifies  the  rare  isotope.  For  example,  in  the  case  of  12C-13C, 
the  uncertainty  in  the  ratio,  R  =  12C/13C,  is  then  given  by  +RAA'13 
where  N 13  is  the  total  number  of  particles  detected  for  the  process 
which  identifies  the  presence  of  C.  For  =  10  the  uncertainty 
in  the  ratio  is  then  +  0.1%.  In  practice  the  total  number  of  detected 
particles  will  be  limited  by  the  cross  section,  maximum  beam  current 
which  the  target  can  withstand,  the  magnitude  of  the  background,  and 
the  total  time  available  for  the  measurement. 

III.  Determination  of  R  by  Elastic  Scattering 

As  we  have  seen  in  equation  2,  the  elastic  scattering  process  ca 
be  used  to  differentiate  between  different  target  nuclei  by  a  simple 


measurement  of  the  energy  of  the  outgoing  particle.  The  mass  discrimina¬ 
tion  is  best  for  180°.  Although  the  scattering  cross  sections  will 
be  largest  at  forward  angles  this  advantage  is  negated  by  the  smaller 
energy  separations  of  the  scattered  particles.  A  rough  upper  limit 
to  target  thickness,  t,  that  may  be  used  without  seriously  degrading 
the  system  can  be  estimated  by  letting  the  maximum  energy  loss  in 
the  target,  dE/dx  •  2t  for  back  angle  scattering  be  equal  to  or  of 
the  order  of  the  energy  resolution  of  the  detector  used.  Since  scattering 
can  occur  at  the  front  or  rear  surface  of  the  target  this  energy  loss  is 

the  energy  spread  of  the  scattered  beam.  For  4  MeV  protons  in  polyethylen 

2 

for  example,  dE/dx  is  100  keV/mg/cm  and  the  resolution  of  commonly  used 

silicon  detectors  is  ~  20  keV;  thus  the  maximum  organic  target  thickness 
2 

is  '^0.2  mg/cm  . 

The  scattering  ratio  for  protons  scattered  by  targets  with  atomic 

masses  from  12  to  18  is  shown  in  Figure  1.  Figure  2  shows  the  range  of 
■ 

values  obtained  for  incident  particles  with  atomic  masses  of  1,  2,  3, 

•  and  4 . 

Now, from  reference  to  Table  I,  it  can  be  seen  that  the  nuclear 

reactions  produced  by  beams  of  protons  and  alphas  are  minimal  compared 

3 

to  those  produced  by  deuterons,  tritons,  and  He.  The  latter  particles 
produce  many  unwanted  reaction  particles  which  complicate  the  detection 
of  the  elastically  scattered  particles  and  hence  are  not  as  useful  as 
protons  and  alphas.  It  is  not  so  easy  to  choose  between  protons  and 
alphas,  however.  Alphas  produce  a  bigger  kinematic  shift,  but  dissipate 
more  energy  in  the  target  and  are  somewhat  harder  to  detect  with  good 
resolution  in  semiconductor  particle  detectors* 


For  a  particular  example  we  will  take  the  case  of  the  measurement 

12  13 

of  C-  C  ratios.  A  sample  with  the  composition  of  polyethylene  and 

a  proton  beam  current  of  1.0  ^A  are  assumed.  The  first  approach  might 

be  to  use  a  reasonably  high  bombarding  energy  ~3.0  MeV  for  example,  in 

order  to  get  as  large  as  possible  difference  in  energies  between  the 

particles  scattered  from  the  two  isotopes.  Target  thickness  has  been 

adjusted  so  that  the  total  rate  is  50,000  counts/sec  which  is  a  reasonable 

limit  imposed- by  commonly  available  electronics.  An  annular  detector 

with  an  acceptance  angle  of  +1  was  assumed  with  an  angle  of  observation 

of  160  .  The  energy  separation  between  protons  scattered  from  the  two 

isotopes  is  then  55  keV .  The  experimental  conditions  with  a  target 

2 

thickness  of  75  ^g/cm  are  such  that  it  should  be  feasible  to  resolve 

13 

the  peaks  corresponding  to  the  two  isotopes.  A  rate  for  c  of  500  to 
1000  counts/sec  can  be  achieved  and  hence  the  abundance  ratio,  R,  can 
be  determined  to  +17.  in  10  to  20  sec  and  to  +.17.  in  1000  to  2000  sec. 

12 

The  large  difference  in  the  counting  rates  for  scattering  from  C 

13 

and  C  also  illustrates  a  possible  limitation;  namely,  the  counting 
equipment  can  be  rate  limited  by  particles  scattered  by  the  abundant 
isotope,  In  order  to  reduce  this  rate  it  is  possible  only  to  search 
for  an  energy  where  the  cross  section  for  the  abundant  isotope  is  a 
minimum  relative  to  the  other  cross  section.  Some  advantage  can  be 
had  from  this,  i.e.  where  non-Rutherford  nuclear  effects  are  operative 

at  low  energies,  but  it  is  not  generally  applicable.  However,  for  the 

12  13  12  12 

C-  C  case  there  is  a  resonance  in  C(p,p)  C  scattering  at  E  ~  1 . 7  MeV 

and  a  rather  narrow  minimum  occurs  there.  If  an  estimate  is  done  for 

13 

this  energy  the  rate  for  C  scattering  can  be  increased  from  the  above 
by  a  factor  of  3  to  5  and  R  can  be  measured  to  +0.17.  in  a  time  of 
perhaps  500  sec. 


For  alpha  particle  scattering  by  C  examination  of  the  excitation 
curves  shows  a  deep  minimum  at  about  E^,  =  3.2  MeV  and  an  angle  of  ~105° 

in  the  laboratory  system.  Choosing  this  energy  and  a  target  thickness 

2  +  • 
of  10  p.g/cm  and  again  a  1  pA  current  of  He  gives  a  rate  of  470/sec. 

12  13 

Here  the  dynamic  range  for  scattering  from  C  and  C  is  very  good;  but 
it  is  probably  difficult  to  increase  the  rate  very  much  by  increasing 
the  beam  current  because  of  the  increased  specific  ionization  and  energy 
loss  to  the  target.  It  again  appears  that  a  determination  of  R  to  +0.1% 
in~2000  sec  is  feasible. 

It  is  concluded  that  it  should  be  feasible  to  measure  the  isotopic 
12  13 

abundance  of  C-  C  to~j^3.1%  in  a  time  of  «%2000  sec  by  elastic  scattering 

techniques.  Better  precision  is  possible  with  longer  counting  periods, 

1/2  ; 

but.it  improves  only  as  t 

14  15 

Similar  considerations  will  apply  to  the  determination  of  N-  H  • 

16  17  18 

and  0-  0-  0  abundances.  In  these  cases  the  overall  composition  of 

the  target  will  be  important  since  the  total  counting  rates  will  probably 
be  determined  by  scattering  from  carbon.  If  nitrogen  and  oxygen  are 
present  in  small  amounts  the  determination  will  be  relatively  inaccurate, 
perhaps  +17.  in  the  ratio,  for  counting  periods  of  2000  sec.  If  a  more 
sophisticated  magnetic  spectrometer  was  used  it  would  be  possible  to-' 
eliminate  the  unwanted  particles,  but  on  the  other  hand  it  will  be  at 
the  expense  of  solid  angle  so  that  the  final  result  may  be  similar. 


IV.  Determination  of  R  by  Nuclear  Reactions 

The  advantage  of  determining  R  by  the  use  of  a  suitable  nuclear 


reaction  is  that  one  need  not  be  limited  by  the  rate  from  the  abundant 


isotope,  but  can  choose  a  reaction  which  will  optimize  the  detection  of 

the  rare  isotope.  While  m'any  of  the  charged  particle  reactions  give 

reaction  products  which  can  be  identified  uniquely  the  cross  sections 

are  in  general  quite  low.  Neutron  producing  reactions  in  several 

cases  present  very  attractive  features.  The  cross  sections  can  be  quite 

large  and  are  often  substantially  bigger  than  the  cross  sections  for 

reactions  producing  charged  particles.  The  neutrons  can,  by  choice 

of  suitable  detection  systems,  be  detected  with  high  efficiency  large 

solid  angle  detectors.  The  neutrons  produced  by  reactions  in  other 

cons itituents  of  the  target  can  again  be  rejected  by  proper  choice  of 
type  of  detector  and  bombarding  energy.  For  these  reasons  only  reactions 

which  produce  neutrons  will  be  considered. 

The  neutron  producing  reaction  is  used  to  indicate  the  abundance 
of  the  rare  isotope  and  the  elastic  scattering  technique  can  be  used 
to  give  either  the  amount  of  the  abundant  isotope  or  the  total  amount 
of  the  element  present.  Equation  5  must  now  be  multiplied  by  the  ratio 
of  the  solid  angle-efficiency  product  for  the  two  detectors.  For  a 
particular  experimental  arrangement  relative  values  of  R  can  be  measured 
accurately,  but  absolute  values  will  be  rather  imprecise.  This  emphasizes 
the  need  to  use  a  standard  sample  for  calibration. 

For  the  12C-13C  measurement,  inspection  of  Table  I  shows  that  the 
13C(a,n)160  reaction  is  attractive  since  neutrons  produced  by  the  alpha 
particle  bombardment  of  the  other  light  nuclei  can  be  easily  discriminated 
against  with  a  simple  hydrogeneous  proton  recoil  type  neutron  detector. 
Making  the  rather  conservative  assumption  that  the  neutrons  are  detected 
with  57.  efficiency  in  a  solid  angle  of  .09  sr.  and  that  a  100  kcV  thick 


target  can  be  used  at  an  alpha  energy  of  3.65  MeV  where  the  cross 
section  is  30  mb/sr,  then  10*  neutrons  can  be  detected  in  a  time 
of  about  300  sec.  The  time  required  to  measure  R  to  ±  0.17.  is  thus 
less  than  the  time  for  the  same  determination  by  elastic  scattering 
by  a  factor  of  about  5. 

For  the  1AN-15N  pair  a  similar  technique  can  be  used  with  the 
(d,n)  reaction.  'The  cross  section  is  not  very  large,  about  5  mb/sr 
at  2  MeV  which,  however,  can  be  compensated  for  by  a  thicker  target. 

The' use  of  a  thicker  target  would  spoil  the  energy  resolution,  but 
the.UN(d,Po)15N  reaction  rather  than  elastic  scattering  could  be  used 
to  measure  the  amount  of  present.  Since  the  (d,pQ)  reaction 

a  O-value  much  more  positive  than  the  other  possible  deuteron  induced 
charged  particle  reactions  no  great  resolution  in  the  charged  particle 
detector  is  needed  and  adequate  rates  can  be  obtained.  The  counting 
time  required  to  make  a  determination  of  the  ratio  would  depend  on 
the  amount  of  nitrogen  present  in  the  sample  of  interest. 

The  case  of  the  160-180  pair  can  be  handled  by  use  of  the  (p,n) 

reaction.  180  is  the  only  light  isotope  which  produces  neutrons  if  the 

proton  energy  is  kept  below  3  MeV.  A  target  thickness  and  angle  of 

observation  are  chosen  which  permit  a  semiconductor  detector  to  resolve 

protons  elastically  scattered  by  160  from  protons  scattered  by  other 

elements.  The  determination  can  then  be  made  at  a  proton  energy  of 

2.650  MeV  where  there  is  a  maximum  in  the  0(p,n)  F  cross  section. 

The  average  total  cross  section  is  25  mb.(11)  The  neutrons  have  low 

energy  and  can  be  detected  in  a  properly  designed  neutron  counter  with 

* 

essentially  1007,  efficiency  and  a  solid  angle  close  to  Arr  sr 


If  a 


51 


2  • 
target  of  LiOtl  17  pg/em  thick  Is  assumed,  which  results  m  a  proton 

energy  loss  of  10  keV,  the  resulting  rate  is  8xl05/see  and  the  de¬ 
termination  can  be  made  accurately.  The  actual  counting  rate  would 
be  at  least  a  factor  of  10  less  because  of  rate  problems  in  the  de¬ 
tection  system. 

17 

The  measurement  of  0  by  the  reaction  technique  does  not  look 

as  promising  as  for  the  other  isotopes.  The  (d,t)  reaction  might 

be  of  use  but  would  require  the  use  of  particle  identification  techniques 

which  are  cumbersome.  The  (o'.n)  reaction  is  perhaps  the  most  hopeful,  but 
13  17 

here  since  the  0(0", n)  F  reaction  has. a  more  positive  0-value,  time-of- 

flight  techniques  would  be  needed  for  a  rather  precise  neutron  energy  measurement. 


V. 


Examples  of  160-180  and  l2C-13C  Tsotopic  Abundance  Ratios  Determined 


by  Elastic  Scattering 

In  order  to  illustrate  the  application  of  the  elastic  scattering 

technique  several  spectra  were  obtained  for  the  energy  distribution 

of  protons  elastically  scattered  by  targets  containing  C  and  0.  The 

angle  o  observation  was  about  160°  and  the  counter  solid  angle  was 

1,65x10"**  steradians.  The  equipment  that  was  used  was  designed  for 

other  experiments  and  no  attempt  was  made  to  minimize  data  collection 

times  for  this  particular  use. 

The  first  target  that  was  used  represents  one  possible 

technique  whereby  thin  oxygen  targets  ear.  be  prepared  in  an 'easy  way 

1 G  18 

for  the  determination  of  0-  0  ratios  in  water  samples.  A  thin  Li 

2  2 
target,  100  pg/cm  ,  was  vacuum  evaporated  onto  a  20  pg/cm  carbon 

foil.  The  target  was  then  removed  from’the  evaporator  and  the  lithium 


allowed  to  either  oxidize  or  form  the  hydroxide  by  reacting  with 

\ 

the  oxygen  and  water  vapor  present  in  the  atmosphere.  This  technique 

gave  a  stable  target  with  oxygen  as  the  heaviest  constituent. 

This  method  could  be  used  in  general  to  prepare  thin  oxygen 

targets  from  samples  of  oxygen  gas  or  water.  The  fresh  Li  target 

would  be  exposed  directly  to  the  oxygen  or  over  the  water  sample  of 

interest  without  exposure  to  atmospheric  oxygen.  Bubbling  of  a  noble 

or  non-reacting  carrier  gas  through  the  water  could  hasten  the  reaction. 

The  targets  should  be  totally  converted  from  the  metallic  form  and 

should  be  stable  enough  to  be  transferred  after  preparation  to  the 

scattering  chamber  with  minimal  precautions. 

In  Fig.  3  data  is  shown  for  proton  bombardment  of  a  Li-C-0 

target  with  a  beam  current  of  about  0.5  uA  for  two  hours.  In  this 

13  18 

case  the  statistical  accuracy  for  the  C  peak  is  ~+47..  The  0 

peak  is  not  as  strong  and  the  statistical  accuracy  is  ~+107..  The 

18  13 

determination  of  the  0  abundance  is  more  difficult  than  for  JC 

because  the  relative  isotopic  abundance  is  less. 

Figure  4  shows  two  pulse  height  spectra  obtained  for  proton 

bombardment  of  two  carbon  foils.  The  target  used  for  the  spectrum 

2 

shown  in  Fig.  4a  was  a  20  pg/cm  foil  of  natural  carbon.  For  Fig. 

2  13 

4b  a  15  pg/cm  foil  enriched  in-  C  was  used.  Analysis  of  this  data 

13 

gives  the  composition  of  the  enriched  foil  as  (84  +  8%)  C  relative 
to  an  assumed  natural  value  of  1.117..  It  should  again  be  emphasized 
that  the  data  shown  here  was  taken  to  il lustrate  the  technique  and 
no  effort  was  made  to  minimize  data  collection  times  or  to  attain 
high  precision.  .• 
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VI.  Conclusions 

1 

We  have  examined  the  possibility  of  determining  isotopic  abundances 
of  the  C,  N,  and  0  isotopes  in  organic  samples  largely  free  of  heavy 
elements.  Very  simple  elastic  scattering  measurements  should  be  capable 
of  determining  the  ratios  to  about  +0.1/1  where  the  element  of  interest 
is  a  major  constituent  of  the  sample.  In  some  cases  the  use  of  a  nuclear 
reaction  to  tag  the  rare  isotope  of  interest  may  give  somewhat  improved 
results.  Determination  should  be  possible  with  running  times  of  5-30 
min  so  that  the  analysis  of  many  samples  should  be  convenient. 

This  method,  while  it  may  not  have  the  ultimate  sensitivity  of 
mass  specurometric  methods,  does  have  several  advantages.  It  does  not 
depend  on  the  cleanliness  of  an  ion  source  and  therefore  may  be  more 
suitable  for  processing  large  batches  of  samples.  Possible  preferential 
extraction  of  one  isotope  by  the  ion  source  is  not  a  problem.  Targets 
can  be  prepared  from  small  amounts  of  material,  perhaps  1  to  10  mg. 

It  appears  that  the  use  of  small  Van  de  Graaff  accelerators  for 
routine  analysis  of  C,  N,  0  isotopic  concentrations  shows  promise 

of  convenience,  ease  of  determination  and  fair  accuracy  and  could  be  a 
valuable  tool,  for  isotopic  analysis.  Even  though  some  progress  has 
been  made  in  such  applications „ . it  appears  that  the  field  has  not 
yet  been  well  explored.  Close  cooperation  between  biologists  and 
chemists  and  others  who  are  interested  in  the  results  Of  specific 
measurements  with  the  experimental  nuclear  physicist  who  is  knowledge¬ 
able  about  the  experimental  equipment  will  be  necessary  to  most  fruit¬ 
fully  exploit  these  basically  nuclear  techniques. 


56 

Acknowledgments 

We  wish  to  thank  Drs .  Leona  Libby,  Martin  Posner,  and  Harvey 

\ 

Wegner  for  helpful  discussions  and  suggestions. 


Table  I.  Q-v 

alues  in 

kcV  for  re 

actions  induced  bv  p, 

o  ,  t,  3i 

i!e,  and 

*He 

arc 

bombardment  o£  stable  isotopes 
taken  from  ref.  6. 

of  C,  N, 

0.  The 

values 

12 

action/Target  C 

'l3C 

. 

15» 

■,16° 

17° 

18o 

p,n 

-18390 

-  3004 

-  5931 

-  3543 

-16212' 

-  3544 

-  2450 

p,d 

.  -16496 

-  2722 

-  8329 

-  8610 

-13444 

-  1917 

-  5822 

P,3 

-23330 

-15186 

-22400 

-12905 

-20408 

-11329 

-  3706 

P,,He 

eM 

-19694 

-13185 

-  4779 

-10666 

-15242 

-  8550 

-14097 

-  7557 

-  4064 

-  2922 

4964 

-  5218 

1193 

3980 

d,n 

.  -  281 

5325 

5066 

9901 

-  1627 

3372 

5768 

d,p 

‘  2722 

5951 

8610 

275 

1917 

5822 

1732 

d,3 

-12463 

1311 

-  4296 

-  4577 

-  9411 

2116 

-  1789' 

d ,  He 

-10463 

-12040 

-  2056 

-  4715 

-  6633 

-  8275 

-10570 

d/4He 

-  1342 

5167 

13574 

7687 

3110 

9802 

4256 

t.n 

.  4014 

9901 

14478 

7786 

1257 

7557 

6109 

t,p 

4641 

911 

4852 

-  230 

3706 

3521 

31.33 

M. 

-  1311 

1918 

4  577 

-  3758 

-  2116 

1789 

-  2301 

*•> 

-13351 

-13410 

-  138 

-  9755 

-10390 

-  8780 

1  3984 

t.  He 

3857 

2279 

12263 

9605 

'  7687 

6045 

3750 

3Hc  ,  n 

-  1152 

7123 

969 

5006 

-  2980 

4290 

13118 

3He,p 

4779 

10666 

15242 

8550 

2021 

8321 

6873 

^He,d 

3He,t 

..  -  3550 

2056 

1798 

6633 

-  4895 

104 

2499 

-17620 

-  2239 

-  5166 

-  2779 

-15448 

-  2780 

-  1686 

JHe/JHe 

1856 

15630 

10024 

9743 

4908 

16435 

12530 

-  8508 

2214 

-  4737 

-  6430 

-12145 

588 

-.  700 

>»P 

-  4964 

-  7411 

-  1193 

-  3980 

-  8114 

-  5657 

-  5620 

lle.d 

-13574 

•  7687 

-  3110 

-  9802 

-16332 

-10031 

-11479 

**He,3He 

*He,3-He 

-17869 

-12263 

-12522 

-  7687 

-19215 

-14216 

-11820 

-15630 

-12401 

-  9743 

-18077 

-16435 

-12530 

-16621 

REFERENCES 


1.  See  for  example:  S.  Brenner,  F.  Jacob,  and  H.  Meselson,  Nature  190, 

•  576  (1961);  M.  Meselson  and  J.  J.  Weigle,  Proc.  Nat.  Acad.  Sci.  47 , 

857  (1961);  Carl  L.  Schildkraut,  Julius  Marmur,  and  Paul  Doty, 

J.  Mol.  Biol.  3,  595  (1961);  Matthew  Meselson  and  Franklin  W.  Stahl, 
Proc.  Nat.  Acad.  Sci.  44,  671  (1958);  Mildred  Cohn;  Journ.  Biol. 

Chem.  201,-735  (1953). 

2.  J.  B.  Marion  and  F.  C.  Young,  Nuclear  Reaction  Analysis,  Korth-Holland , 
Amsterdam,  1968. 

3.  G.  Amsel ,  J.  P.  Nadai,  E.  D'Artemare,  D.  David,  E.  Girard,  and 
J.  Moulin,  Nucl.  Inst.  92,  481  (1971). 

4.  B.  L.  Cohen  and  R.  A.  Mozer,  Anal.  Chem.  43,  123  (1970. 

5.  R  K.  Jolly,  C.  R.  Gruhn,  and  C.  Maggiore,  IEEE  Trans.  Nucl.  Sci. 

NS-18 ,  91  (1971). 

6.  F.  Ever ling,  L.  A.  Koenig,  J.  H.  E.  Mattauch,  and  A.  H.  Wapstra,  I960 
Nuclear  Data  Tables,  Part  1,  Nuclear. Data  Project,  NAS-NRC  (1961); 

A.  H.  Wapstra  and  N.  B.  Gove,  Nuclear  Data  Tables  9,  267  (1971). 

7.  Los  Alamos  Report  LA-2014  (1957),  edited  by  Nelson  Jarmie  and 
John  D.  Seagrave. 

8.  H.  J.  Kim,  W.  T.  Milner,  and  F.  K.  McGowan,  Nuclear  Data  A  2,  1  (1966). 

9.  F.  Ajzenberg-Selove  and  T.  Lauritsen,  Nucl.  Phys . ,  A  114,  1  (1968). 

10.  F.  Ajzenberg-Selove,  Nucl.  Phys.  A152,  1  (1970). 

11.  J.  M.  Blair  and  J.  J.  Leigh,  Phys.  Rev.  118,  495  (I960). 


59 


Fig. 
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Fig.  3 


Fig.  4 


FIGURE  CAPTIONS 

Scattering  ratios  for  protons  scattered  by  target  nuclei  with 
masses  between  12  and  18  as  a  function  of  laboratory  scattering 
angle.  Eq  is  the  energy  of  the  incident  particle  and  is  the 
energy  of  the  scattered  particle. 

Scattering  ratios  for  mass  1,  2,  3,  and  4  bombarding  particles 
scattered  by  target  nuclei  with  masses  between  12  and  18  as  a 
function  of  laboratory  scattering  angle.  Target  nuclei  with 
intermediate  masses  have  scattering  ratios  which  fall  in  the 
cross  hatched  areas.  E^  is  the  energy  of  the  incident  particle 
and  E^  is  the  energy  of  the  scattered  particle. 

Portion  of  a  pulse  height  spectrum  showing  the  scattering  of  3.0 
HeV  protons  by  a  target  of  Li,  C,  and  0.  The  laboratory  angle  of 
observation  is  160°. 

The  spectrum  shown  in  4a  shows  the  elastic  scattering  of  3.0 

o  2 

MeV  protons  at  an  angle  of  156  by  a  20  ^g/cm  natural  carbon 

foil.  In  4b  the  spectrum  obtained  under  identical  conditions 
2  13 

with  a  15  pg/cra  C  enriched  foil  is  shown.  The  13C  abundance 
in  the  enriched  foil  was  found  to  be  (84  +8)7,  from  this  data. 


